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Abstract. We present new elements in the identification 
of the lens-system producing the 4 images of the BAL 
quasar H1413+117, based on the recent HST/NICMOS- 
2/F160W observations. After a careful PSF subtraction 
of the 4 images of the quasar, the residual Hfi60W im- 
age reveals the presence of a faint object (H ~ 20.5) 
within the region enclosed by the 4 quasar images. This 
object corresponds to a single galaxy: the primary lens 
of the lens-system. We also identify the galaxies around 
the Cloverleaf which had been proposed to belong to a 
galaxy cluster/group at high redshift (Kneib et al 1998): 
the other component in the lens-system that provides the 
additional "external" shear. For these galaxies, we have 
derived a likely redshift based upon their Kf702W, Ifsi4W 
and Hp i6oiv magnitudes. We find that most of them are 
consistent with belonging to a galaxy cluster/group with 
mean redshift z = 0.9 ±0.1. Furthermore we detect 2 
very red objects (I — H ~ 4): the faintest one has no ob- 
served optical (Rf702W and \fsiaw) counterpart, while 
the brightest has a predicted redshift around z ~ 2, and 
may be identified with one of the Cloverleaf absorbers. 
This gravitational-lens system constitutes an excellent tar- 
get for IR imaging/spectroscopy with the new generation 
of 8m ground-based telescopes. 
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1. Introduction 

The excellent quality and broad wavelength coverage of 
current observations of gravitational lensing systems have 
allowed to unveil part of their mysteries. The primary lens 
is often detected. The immediate surrounding of the multi- 
ply imaged quasars/galaxies is studied in great detail and 
generally shows some galaxy clustering (e.g. Tonry 1998, 
Hjorth & Kneib 1998) or even sometimes X-ray cluster 
emissions (Hattori et al 1997, Chartas et al 1998). Further- 
more, the measure of a time-delay in Q0957+561/PG1115 
(Kundic et al 1997, Schechter et al 1997) has strengthened 
the interest of a detailed study of these gravitational lens 
systems in order to use them as a cosmological tool. 

Since the identification by Magain et al (1988) of the 
quadrupole lens-system called the Cloverleaf, 4 images of 
the BAL quasar H1413+117 at z=2.558, many efforts have 
been dedicated to a direct search of the lens or of elements 
of the lens-system. Early models of the lens-system have 
involved one or two galaxy-lenses very close to the line-of- 
sight toward the quasar (Kayser et al, 1990). A more re- 
cent analysis showed that an external shear was needed to 
model this system correctly (Keeton, Kochaneck & Seljak 
1997), and indeed it is probably related to the existence 
of an overdensity of galaxies nearby, as detected by Kneib 
et al (1998). 

The lensing geometry, amplification and time-delays 
are quite sensitive to the parameters of the galaxy-lens 
expected to be located amid the 4 images of the quasar. A 
positive detection of the galaxy-lens would bring stringent 
constraints in the modeling of the lens-system. However, 
despite relatively deep searches, either in R and I imaging 
with the HST (Turnshek et al, 1997; Kneib et al, 1998) or 
K imaging with the Keck telescope (Lawrence et al 1996), 
the lensing galaxy has not been detected. 

The galaxy cluster recently revealed near the Clover- 
leaf (Kneib et al, 1998) was assumed to be at z ~ 1.7 as 
this corresponds to the mean value of the narrow absorp- 
tion line-systems observed in the quasar spectra (z=1.44, 
1.66, 1.87, 2.07 and 2.09: Turnshek et al, 1988; Magain 
et al, 1988, Monier et al 1998). Combining the IR image 
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Fig. 1. Final NICMOS-2/F160W image of the Clover- 
leaf. Ellipses indicate the objects detected in the 
WFPC2/F814W image (following the numbering of Kneib 
et al 1998). 



with the F1F702W and Ifsiw WFPC-2 images of this sys- 
tem can allow to estimate the likely redshift for the faint 
galaxies surrounding the Cloverleaf. The recently acquired 
HST/NICMOS-2/F160W observation consists of a unique 
dataset to help answer both the existence of the lensing 
galaxy and to constrain the distance of the nearby galax- 
ies. 

The HST/NICMOS-2 data are presented in Sect. 2, 
while in Sect. 3 we discuss the identification of the lens- 
ing galaxy after the PSF subtraction of the 4 quasar im- 
ages. The redshift estimates (derived from photometry) 
of the surrounding galaxies are explained in Sect. 4. The 
discussion and concluding remarks are provided in Sect. 5. 
Throughout this paper we use Ho —50 km/s/Mpc and 

n Q =i. 



2. The HST/NICMOS-2 data 

We retrieved the NICMOS-2/F160W observations of the 
Cloverleaf which do not hold proprietary rights from the 
HST/STSci archive: http://archive.stsci.edu/cgi-bin/hst. 
These data are part of a much larger dataset (PI: E. Falco) 
with the goal of studying the statistical properties of a 
sample of gravitational lens systems. The Cloverleaf data 
were taken on December 28, 1997 with a total exposure 
time of 4x640 seconds. At the date of the observations, 
the NICMOS-2 pixel size (x,y) was 0.07603" x 0.07534" 
(c.f. the NICMOS instrument WEB-page). The NICMOS 



data come in 3 different formats: raw (files *_raw.fits), 
calibrated (files *_cal.fits) and mosaiced (files *_mos.fits). 
The raw data set is made of 4 sets of 19 rasters. The cal- 
ibrated data correspond to 4 flux calibrated images, each 
being a combination of the 19 flat-field corrected images. 
We made our own mosaiced image by combining the 4 
flux calibrated images shifted with a non-integer number 
of pixels (Figure ^) . The image quality measured on the 
final mosaiced image is 0.145" (FWHM). 

The photometric calibration in the Vega system was 
determined from the PHOTNU image keyword following 
the NICMOS handbook recipe (PHOTNU=2.3857 10" 6 ). 
The software SExtractor 1.2 (Bertin & Arnouts 1996) was 
used to detect and to derive the photometry of faint sources, 
the detection parameters used are: 15 contiguous pixels 
over 1.5cr of the sky-background after a convolution with 
a 3x3 top- hat filter. Photometry and position relative to 
the brightest quasar image (A) of the different objects is 
summarized in Table |^. 

3. The 4 quasar images and the primary 
galaxy-lens 

In order to reveal the primary galaxy-lens amid the 4 
quasar images, we have subtracted modeled PSFs of the 
NICMOS-2 camera. We used the TinyTim V4.4 software 
(Krist & Hook 1997) to model the NICMOS-2 PSF through 
the F160W filter (assuming a flat spectrum for the quasar 
through this waveband) at each of the 4 positions of the 
quasar images and in each of the 4 calibrated frames. We 
oversampled the PSF with 8x8 subpixels per NICMOS-2 
pixel (in order to sample correctly the PSFs and position 
them accurately). The positions of the 4 images of the 
Cloverleaf given in Turnshek et al (1997) were used for the 
PSF subtraction (we also independantly fitted the position 
directly from the NICMOS image but it did not improved 
the substraction), while the flux ratios were adjusted in 
order to minimize the mean residuals at the 4 quasar im- 
age locations. This procedure was applied to each of the 
4 calibrated images, allowing us to derive the magnitudes 
and flux ratios of the 4 quasar images in Hfwow (Table 
0), and it also gives independant errors on the photome- 
try of the different images. We found systematic residual 
patterns for each of the 4 PSF-substracted quasar images 
- as it can be seen in Figure 2 - but these patterns do not 
corresponds to the location of the detected lensing galaxy. 
Hence we could not check for the existence of an Einstein 
ring. Flux ratios in H differ significantly from the ones 
in Uf336w, Rf702W and I_F8i4VK, particularly for image 
D. As proposed by Turnshek et al (1997), dust extinction 
can easily explain the different flux ratios. However mi- 
crolensing on the D image is also very likely, as it is the 
closest image to the detected lensing-galaxy. 

In order to quantify the magnitude and position of the 
galaxy-lens, we assume that its profile follows a de Vau- 
couleurs law and we include it in the fitting of the PSF 
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Fig. 2. PSF subtraction steps, from left-to-right and top- 
to-bottom: the NICMOS-2 observed image, schematic of 
quasar-images and galaxy positions, PSF model for the 
4 quasars, PSF subtraction of the 4 quasars, PSF model 
for the 4 quasars plus the galaxy-lens, PSF subtraction of 
the 4 quasars and the galaxy- lens. Size of these rasters is 
3.4 arcseconds on a side. The orientation corresponds to 
the NICMOS-2 images (PA of image from Y axis 74.3765 
degree) . 



subtraction. Parameters were chosen in order to be left 
with a mean sky value similar to the median sky value. 
However we could not constrain the structural parame- 
ters of the galaxy-lens due to poor signal to noise. Figure 
^ exhibits the different steps in the subtraction procedure. 
Note that the galaxy-lens, HI, is clearly detected even on 
the non-PSF subtracted image. Its position and Hfi60W 
magnitude are given in Table |^. Unfortunately, we can 
only get lower limits to its Rf702W and Ifsuw magni- 
tudes, preventing us from deriving a redshift constraint 
for HI. 



4. The environment of the Cloverleaf: 
photometric redshifts of the faint objects in the 
field 



In the field-of-view covered by the Hfwow image, we iden- 
tify 12 faint objects for which we can provide at least one 
of the Rf702W, I.F814W or Hfwow magnitude [in addition 
to HI]. Among these objects, 11 have been identified by 
Kneib et al (1998). Another one, called H2, was not de- 
tected either on the Rf702W or Ifsuw images. All these 
objects are shown on Figure [|. Their positions with re- 
spect to the quasar image A are given in Table ||, as well 
as their Rf702W: Ifsmw and Hfwow magnitudes. 

We have estimated the redshift of these objects from 
their Rf702W, Ifsmw and Hfi60W magnitudes, accord- 
ing to the standard minimization method described by 
Miralles & Pello (1998). The updated Bruzual & Chariot 
evolutionary code (Bruzual & Chariot, 1993) was used and 
5 different synthetic star-formation rates (SFR) tested: 
a 0.1 Gyr burst, a constant SFR, and 3 \i models (e- 
decaying SFR) with characteristic decay-time matching 
the sequence of colors for E, Sa and Sc galaxies. For each 
SFR type, 51 template spectra were selected with differ- 
ent stellar population ages. We also tested the stability of 
our results against different Ay values in the galaxies. The 
template database includes 255 solar metallicity spectra. 

We have analyzed through simulations the accuracy 
and the possible biases affecting the redshift estimate from 
the photometry. For this purpose, a simulated catalogue of 
galaxies was created, with redshifts uniformly distributed 
between and 6, and randomly sampling the SFR and 
ages. Photometric errors ( uncorrelated for the different fil- 
ters) were introduced as gaussian noise distributions with 
FWHM chosen so as to match the typical values for the 
Cloverleaf galaxy sample. According to these simulations, 
the use of the Rf702W, Ifsmw and Hpieow filters does 
not introduce any systematic bias in the redshift determi- 
nation for 0.5 < z < 2.5. 

The most probable redshift was obtained for each ob- 
ject, as well as its probability (in the sense of x 2 ) as a 
function of redshift (Table |2|). We find that 9 objects have 
a redshift between 0.8 and 1.1, with z = 0.9±0.1. Most of 
these objects are well fitted by a 1 to 2 Gyr burst models 
and Ay — 0. The If814W filter maps the B restframe at 
z ~ 0.9. The absolute magnitudes of these galaxies are 
-2K M B <-19. 

For 4 objects (including the galaxy-lens HI and the 
very red object H2) the redshift is poorly determined. 
The objects #14 and H2 are extremely red in (If814W _ 
Hfi6(w) and are consistent with galaxies at redshift 2 
or higher. For object #14, the best-fit model (without 
reddening), is a 2.5 Gyr burst, which, given the redshift 
uncertainty, leads to Mb = —22.5 to —23.9. This high lu- 
minosity may suggest that object #14 might as well be 
gravitationally amplified by the cluster- lens with z ~ 1. 
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Uf336VK 










(") 


K > 


94/12/23 


94/12/23 


94/12/23 


97/12/28 


A 


0.000 


0.000 


18.95 


17.68 


17.45 


15.67 


B 


0.744 


0.172 


19.31 [0.73] 


17.82[0.88] 


17.58[0.89] 


15.76[0.92] 


C 


-0.491 


0.716 


19.93 [1.01] 


17.95[0.77] 


17.73[0.76] 


14.03[0.71] 


D 


0.355 


1.043 


19.13 [0.85] 


17.97[0.74] 


17.81[0.70] 


14.29[0.57] 



Table 1. Relative positions (from Turnshek et al 1997), photometry (Uf336W, Rf702W A-fsuw & Hi?i6oiy) and flux 
ratios of the 4 quasar images relative to image A. Note the strong variation of the flux-ratio from U to H. 
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max 


75% 




HI 


0.12 


0.50 


>22.5 


R >22.7 


20.5±0.7 






the position accuracy is ±0.04 arcsec 


H2 


-7.6 


-7.2 


>24.8 


R >25.5 


21.6±0.7 






very red, only detected in H. 


1 


-7.3 


-2.1 


22.52 


0.72 


20.7±0.3 


0.8 


0.6-1.1 


E/S0 morphology 


2 


-5.5 


-3.7 


23.11 


0.86 


22.3±0.9 


0.9 


0.8-1.0 




3 


-0.6 


8.2 


23.16 


0.88 


22.0±0.5 


0.9 


0.8-1.0 


disturbed morphology 


7 


3.9 


-5.5 


23.75 


0.74 


>22.3 


0.9 


0.5-1.2 




8 


-0.7 


-4.1 


23.78 


1.16 


21.8±0.5 


1.0 




z p hot solution at ~ 50% confidence level 


9 


-5.2 


2.2 


23.94 


0.12 


>22.3 






no Zphot constraint 


11 


3.2 


-6.8 


24.01 


1.09 


22.3±0.8 


1.0 


0.7-1.2 




14 


4.7 


0.2 


24.37 


>1.1 


20.0±0.2 


2.0 


1.8-2.7 


very red, possibly associated with absorbers 


15 


-2.2 


-7.8 


24.41 


1.05 


>22.3 


0.8 


0.6-1.2 




19 


-0.9 


5.3 


24.74 


0.88 


>22.3 


1.0 


0.5-1.5 




20 


-6.1 


-6.1 


24.78 


>0.7 


>22.3 


1.1 


0.2-1.6 


Zphot is poorly determined 



Table 2. Relative position (to quasar image A), photometry (Rf702VK J-fsiaw & ^Fimw), and predicted redshift of 
the galaxies within 10 arcsec around the Cloverleaf. The best matching redshift is given, as well as the permitted 
redshift interval corresponding to a 75% confidence level. 



The photometric redshifts were computed using a solar 
metallicity. This is a fair approximation as, for 1< z <2 
and the filters considered in this analysis, the results are 
rather unsensitive to metallicity. Furthermore, the esti- 
mated redshifts given in Table |^ are insensitive to a red- 
dening Ay < 2. If the reddening is much higher, the es- 
timated redshift would shift toward a lower value. How- 
ever there is no reason to believe that all objects in the 
Cloverleaf field should be highly reddened. The dispersion 
is high, however, as reflected in the large interval at the 
75% confidence level (Table 2); it is therefore important 
to effectively measure the redshift of the brightest objects 
to validate this conclusion. 

5. Discussion and concluding remarks 

Two main results have been obtained from the NICMOS- 
2 data. For the first time the galaxy-lens, HI, close to the 
line of sight toward the quasar, has been identified. Its 
position with respect to the quasar line-of-sight is found 
to be similar (within the uncertainties) to the one derived 
in the various Cloverleaf gravitational lens models (e.g. 
Kneib et al, 1998). It remains difficult to estimate the 
redshift of the galaxy-lens HI because the PSF subtraction 
leaves an increased background noise in the region amid 
the 4 quasar images. Yet, a redshift estimate around 1.0 
or higher is consistent with the Hfigow magnitude and 



the Ifsuw lower limit magnitude we have derived for HI. 
Clearly deep spectroscopic data are needed to solve for the 
determination of its redshift. We find also that there is a 
unique galaxy-lens, in contradiction to some early models 
in which two galaxy-lenses had been envisaged (Kayser et 
al, 1990). Assuming that HI is around z ~ 1, and has 
similar colors and absolute magnitude than the nearby 
galaxies gives a Mass-to-light ratio of M(< 5.1kpc)/Ls ~ 
25 M/L B0 . 

With regard to the Cloverleaf environment, we show 
that 8 nearby galaxies have a most probable redshift around 
0.9, giving credit to the presence of a galaxy cluster/group 
along the line of sight to the Cloverleaf. In our previous 
modelling (Kneib et al, 1998), we assumed for this galaxy 
cluster/group a redshift of 1.7, as a mean of the redshifts of 
the 4 absorbers silhouetted on the quasar spectrum. This 
value should be revised. The location of the galaxy-lens is 
now known from the NICMOS-2 observations and will be 
implemented in a new model of the lens-system. One of the 
faint galaxies surrounding the Cloverleaf appears to be at 
a larger redshift, around 2, and might be related with the 
absorber at z=2.07 or 2.09 (Monier et al 1998). Further 
IR imaging/spectroscopy of these galaxies should remove 
the remaining uncertainties of the Cloverleaf lens-system. 
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